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Abstract 

The differential constraints are applied to obtain explicit solutions of nonlin- 
ear diffusion equations. Certain linear determining equations with parameters 
are used to find such differential constraints. They generalize the determining 
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equations used in the search for classical Lie symmetries. 
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q ■ 1. Introduction. 

Differential constraints arisen originally in the theory of partial differential 
equations of the first order. In particular Jacobi used differential constraints to 
find the total integral of nonlinear equation 



F(x\, x n , z, z Xl , z Xn ) 0, 



Konig applied them to the equation of the second order . They required that 
^ ■ the corresponding over-determined system was compatible. The general theory 
of over-determined systems was developed by Delassus, Riquier, Cartan, Ritt, 
Kuranishi, Spencer and others. One can find references in the book of Pom- 
maret Q2J . Now the applications of over-determined systems include such diverse 
fields as differential geometry, continuum mechanics and nonlinear optics. 

The general formulation of the method of differential constraints requires 
that the original system of partial differential equations 

F 1 = 0,..., F m = (1) 

be enlarged by appending additional differential equations (differential con- 
straints) 

/ii = 0,..., ^ = 0, (2) 
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such that the over-determined system (1), (2) satisfies some conditions of com- 
patibility. 

One can derive many exact solutions of partial differential equation by means 
of differential constraints. It was particularly shown in [j3j that some soliton so- 
lutions can be found using differential constraints. Olver and Rosenau Olver 
[|5f, Kaptsov ||, Levi and Winternitz j5f show that many reduction methods 
such as non-classical symmetry groups, partial invariance, separation of vari- 
ables, the Clarkson-Kruskal direct method can be included into the method 
of differential constraints. In practice, methods based on the Riquier-Ritt the- 
ory of over-determined systems of partial differential equations may be difficult. 
The problem of finding all differential constraints compatible with certain equa- 
tions can be more complicated than the investigation of the original equations. 

It was recently proposed a new method for finding differential constraints 
which uses linear determining equations. These equations are more general 
than the classical determining equations for Lie generators |7j] and depend on 
some parameters. Given an evolution equation 

u t = F(t, x, u,ui,..., u n ) , (3) 

where u k = §^|, then according to [Q the linear determining equation corre- 
sponding to (3) is of the form 

N i 

Dt(h) = E E b lk D^\F UN _ k )D*-\h), b lk e R. (4) 

i=0k=0 

Here and throughout D t , D x are the operators of total differentiation with re- 
spect to t and x. Equality (4) must hold for all solutions of (3). The function h 
may depends on t, The number p is called the order of the solu- 

tion of equation (4). If we have some solution h then corresponding differential 
constraint is 

h = 0. (5) 

It was also shown in [j8f that equations (4) and (5) constitute the compatible 
system. Thus we sketch the derivation of some solutions to the evolution equa- 
tion (4). 

(I) Find solutions of the linear determining equations (4). 

(II) Fixing the function h, we obtain differential constrain (5). 

(III) Find the general solution of (5) which includes some arbitrary functions 
di depending on t. 
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(IV) Substitute the general solution into (4). It leads to ordinary differential 
equations for functions at. 

(V) Solve the ordinary differential equations and obtain a solution of the evo- 
lution equation (4). 

In this paper we start with determination of the solutions of linear deter- 
mining equations of the second and third orders for the nonlinear diffusion 
equation 

u t = {u k u x ) x + f(u). (6) 

These solutions exist only if / belongs to the special forms. Then we use the 
obtained functions h to find solutions of equation (6). In final section we derive 
exact solutions of two-dimension equation 

u t = A ln(u). 



2. Solutions of linear determining equations. 

The nonlinear diffusion equation 

u t = (Q(u)u x ) x + f(u). (7) 

often arises in the description of various physical processes. The group classifi- 
cation of the equation has been carried out in [|9fl. Some exact solutions of (7) 
can be found in |TD], [TT|. In physical applications Q is usually taken to be a 
power function. In this section we consider the equation 

u t = (u q u x ) x + f(u), (8) 

where / is an differentiate function, q ^ 0. If q = —2, / = u or / = const 
then the equation (8) can be linearized. We shall not discuss this case here. 
The linear determining equation, which corresponds to (8), is 

D t h = u q D 2 x h + b iq u x u q - l D x h + (hqu q - l u xx + b 2 q{q - l)u q - 2 u 2 x + b 4 f u )h, (9) 

where &i, . . . , 64 G R. We shall seek solutions to (9) in the form 

h = u n + g(t,x,u, . . .,Un-i), 

where n > 2, Uk = ^x- The method for finding solutions is very similar to the 
standard procedure applied in the group analysis of differential equations [|T2 
and only one of all possibilities is described here for the sake of brevity. 
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We set n = 2. First, let us express all ^-derivatives in (9) using (8). As a 
result, the left-hand side of (3.4) becomes a polynomial with respect to 113,112. 
The polynomial must identically vanish. Collecting similar terms we obtain the 
following relations for the coefficients of 113 and u\ 

q{bi - 4) = 0, ug UlUl + q(b 3 - 3) = 0. 

Thus 61 = 4 and g can be represented as follows 

(3 — b 3 )q o , . . . 

g = — — u x + a(u,t,x)u 1 + g 1 {u,t,x); 

here a and g\ must be functions of u, t and x alone. Collecting the coefficients 
of U211I and U2U1, we have the equations 

2b 2 q - 2b 2 - b\q + b 3 q + 46 3 - 6q = 0, (10) 

2ua u + g(6 3 + l)a = 0. 
In follows from the last equation that 

, . Q±M „ 
a = ai(r, x)u 2 % 

where ai is a function of t and x. Next we consider the coefficient u\ and obtain 
equation 

4b 2 q - 46 2 + bjq - 4b 3 q + 26 3 - 9q + 6 = 0. (11) 
From (10) and (11) it follows that 63 = 1 or 63 = 2±2. 

Assuming 63 = 1, we obtain 62 = ^rp The coefficient of U2 give us equation 

u* +2 (2a lx + / tt (6 4 - 1)) + u^ +1 q 9l = 
The equation enables us to express 

91 = -u l - q {f u {l - 64) - 2aiJ 

The coefficient of yields Euler equation 

w 3 (l - b 4 )f uuu + w 2 (2 - qb 4 - 2b A )f uu - uq 2 f u + q 2 f = 0. 

Consider for simplicity the case 64 = 1. It is easy to see that the last equation 
has two types of solutions: 

/ = ku + m/T' 7 , q 7^ — 1 
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or 

f = ku + nuhiu, q = — 1 

where k, n are arbitrary constants. Let us focus on / = ku + nu~ q . It follows 
from above calculations and equation (9) that 

u q+1 (-qa u - 3u q qai xx - ^u q a lxx + kq 2 ai)ui+ 

+nq 2 aiu\ + 2u q+2 (a\ tx - u q a\ xxx - kqa Xx ) + 2una\ x = 0. (12) 

From (12) we have na\ = 0. 

If a\ = then the solution of (9) is h = u 2 + qu 2 /u. If n = and q ^ —4/3 
then one easily computes 



cl\ = (rx + s) exp(kqt) 
,2 , 2 

In the case g = — 4/3, we obtain 
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u 2 ( 2 \ 

h = U2 + q h (rx + s)wi ur )u~ q exp(kqt). 

u \ q J 



a\ = (rx 2 + sx + p) exp(-^-kt), 

o 

. 4u? / 2 3 ,\ 4/ o , 4 

h = Ui — h I (rx + sx + p)ii\ + —u(2rx -\- s) )u exp(--kt). 

Oil \ z / o 

It can be shown that the found functions /i lead to invariant solutions of the 
corresponding equation (8). 

We omit here for the sake of brevity intermediate calculations and give the 
list of solutions to the equation (9): 

(1) if q = — 1 and f = su + ru\n(u) then 

h = u 2 ; 

u 

(2) if q — 1 and f = su + ru~ q then 

7 qui 

h = W2 H ; 

u 

(3) if g = — 2 and / = sit + rw 3 then 

3w? 



h = u 2 - 



2u 
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(4) if q = 1 and / = ru then 

h = ii2 + s exp(rt)w _2 wi + r/3; 

(5) if q is an arbitrary constant and f = su + ru x ~ q then 

(g - l)u\ 

h = u 2 , 

u 

with r, s E R. 

We did not include functions h that correspond to invariant solutions of the 
equation (8). 

If we will look for solutions to the equation (9), which depend on third 
derivative, then obtain the following list: 

(1) if q is an arbitrary constant and / = su + ru l ~ q + n ^ q ^ u q+1 then 



Q 

3(q — 1) . o . u? 

/i = m 3 H wiw 2 + (a - 3g + 2)— + nwi; 

(2a) if g 7^ 1 and f = nu + ^u q+1 then 

(3g — 1) . . Ui ,., . 

/i = u 3 H uiu 2 + q(q - 2)— + rui; (13) 

(2b) if g = — 2 or g = —4/3 and f = nu + L u q+l + mu q+3 then ft is also given 
by (13); 

(3) if g = — | and / = mu then 

h = u 3 - _|_ ^1 _|_ rexp(-3m£/2)w 5/2 + s exp(ra£/2)w 1/2 ; 

(4) if g = — | and f = nu + mu 5 ^ 2 then 

/i = m 3 jt"^ + -T^T + r exp -3nt/2 w 5/2 ; 

2u 4u z 

(5) if g = — \ and / = mu — 2ku l l 2 then 

5miM2 5m? 1/2 
2m — 4^2 + 1 + sex P( m V 2 ) M ; 

(6) if g = — § and f = nu then 

ft = m 3 - + + s exp(-7n/2 V /2 + r exp(-3n*/2) W 5 / 2 ; 



(7) if q = — 1 and / = mu then 

4lliU2 3m? 2 

h = us 1 7T + sexpf — 2mt)u u\\ 

u u z 

with r, s, m, n G i?. Here we also did not include functions h leading to invariant 
solutions of (8). 



3. Solutions of diffusion equations. 

In this section we shall use the functions obtained above to construct so- 
lutions of diffusion equations (8). One can apply the method described in 
introduction. 

We first take the function h = + qu\/u, where q G R, corresponding 
to some cases mentioned above. Simply by equating this function to zero, we 
obtain the differential constraint 

u 2 + qu\/u = 0. (14) 

The equation (14) has two types of solutions: 

u = {c\x + c 2 )9+i, (15) 

u = c\ exp(c2x), q = — 1, (16) 

where ci,C2 are functions of t. 

If we substitute the representation (16) into equation 

Ut = (u x /u) x + ku In u (17) 

then this leads us to differential equations for ci,C2- From this equations it is 
easy to find c\,C2 and obtain the following solution of (17) 

u = s\ exp(s 2 £ exp(Atf)), Si,s 2 G-R. 

Substituting (15) into equation 

u t = (u q u x ) x + su + rz/T' 7 , 

we find the solution 

ax + b — ^ exp(— s(g + l)t)J , 



with a, b G R. 

It is easy to see that the differential constraint 

u 2 + (g- 1)^ = 
u 

for equation 

u t = {u q u x ) x + ru + sw 1-9 

leads to solution 

2 i 

u = (gax exp(gr£) H exp(2gr£) I , a G R. 

Now let us consider the differential constraints of the third order. We start 
with the equation 

m = (u q u x ) x + su + ru 1 ^ + n?-^—u q+1 , neR (18) 

Q 

As explained above this equation is compatible with the differential constraint 

^3 + 3(g - 1)^ + (g 2 - 3g + 2)4 + nu x = 0. (19) 
u u z 

By a change of variable v = u q one may rewrite (18), (19) in the following way 

Vt = vv xx + -vl + n - v 2 + sqv + rg, (20) 

v 3 + nv 1 = Q. (21) 
If n = — 1 then it follows from (21) that 

v(t, x) = a + 6exp(x) + cexp(— x), (22) 

where a, 6 and c are some functions of t. Substituting this representation into 
equation (20) we obtain the system of ordinary differential equations for the 
function a, b and c: 

a t = -a 2 (l + l/g) + asg-46c/g + rg, (23) 

b t = -ab(l + 2/g) + 6sg, (24) 
c t = —ac(l + 2/g) + esq. (25) 
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Using finite-dimensional invariant subspaces, Galaktionov [[T3j found represen- 
tation (22). 

From (24) and (25) we derive the first integral b = kc, k G R. Therefore the 
system (23)- (25) can be reduced to nonlinear ordinary differential equation for 
the function a. In general we can not express solutions of (23)- (25) in terms of 
the elementary functions. We give one example of exact solution of equation 
(20), with q = — 1 and r = 0. This solution has the representation (22) and the 
functions a, 6, c are 

a = cli sin(pri exp(rit) + m) exp(rit) / cos(]9ri exp(rit) + m), 

b = a<i exp(ri^)/ cos(]9ri exp(r\t) +ra), 

c = (Z3 exp(ri^)/ cos(j9ri exp(ri^) +m), 

where a\ = pr\ , 03 = p 2 rf/4a2, ri = — s and p,a,2,m are arbitrary constants. 
Now let us consider the equation 

u t = (u~ l l 2 u x ) x + mu — 2k^/u, m,k G R (26) 

and the differential constraint 



u 3 - ^ + 54 + km + S e^ 2 v^ = 0. (27) 
2w 4ir 



Using the equation (26), one can write (27) as 

(\nu) tx + se mtl2 u- 1 = 0. 

Replacing h\(e mt l 2 ir l ) by w : the last equation is replaced by the Liouville 
equation 

w tx = sexp(w). 
Since the general solution of the Liouville equation is 

2T'X' 

it gives the representation 



9 



where T and X are the arbitrary functions of t and x respectively. Substituting 
this representation into (26), we have 

^1 _ ( T + x)(r')- 3/2 T" - y (T + X){T')- 1/2) ) = 



= -{X')- 3/2 {X") 2 - (X')- l l 2 X"' - 2k{X'fl 2 . 



(29) 



Differentiating (29) with respect to t, it is easy to obtain the equation for T 



2T"'T' -3(T") 2 + ^(T') 2 = 0. 
If m y£ then the function 

(mt \ 

T = ci tanhl — + c 3 I + c 2 , Ci,c 2 ,c 3 Git 

is the general solution of this equation. Substituting the function T into (29). 
we get the following equation for X 



sm ( 



(c x -c 2 -x\= ^{X')- 3/2 {X") 2 - (X'y l l 2 X'" - 2k(X') 1 ' 2 . 



It should be noted that (28) is equivalent to the representation 

u= (ai+a 2 e mi/2 ) 2 , 

where a±, a 2 are functions of x. This representation yields the following system 
for a\ and a 2 

a\ xx + a 2 m/2 - a\k = 0, (30) 

a<2 xx + aia 2 m/2 - a 2 /c = 0. (31) 

In general, it is possible to express a\ in terms of the Weierstrass function p 
and a 2 in terms of Lame's function jp| . However, if m = 12 and k = 4 , then 
the functions 

1 



cosh [xY 



a 



+ 



cosh (x) cosh (x y 
satisfy the equations (30), (31). 



'sinh4x sinh2x 3x\ , „ 

H ^ 1 — — , a,b e R 



32 
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According to our results in the previous section, as k = , the equation (26) 
is compatible with the differential constraint 

u 3 - ^ + ^4 + r e - 3m< /V/ 2 + se mt ' 2 ^ = 0. (32) 

2u 4:U Z 

Using the equation (26), one can write (32) as 

(In u)tx + re- 3mt/2 u + se m */ V 1 = 0. (33) 

Replacing ln(e~ 3m */ 2 w) by w in (33) yields 

w tx + e w + sre- w - mt = 0. 

If we set s = then from the last equation we find the following representation 

2X'T' 

(iTt) 2 ' 



where T and X are the arbitrary functions of t and x respectively. Substituting 
this representation into (26) leads to equation 

^2~fie* mt l*(-rn(T') l ! 2 (X + T) - 2(T / )" 1/2 T // (X + T) + 4(T') 3/2 ) = 

= -2X /// (X / )" 3/2 (X+T) 2 +8X // (X / )" 1/2 (X+r)-8(X / ) 3/2 +(X // ) 2 (X / )" 5/2 (X+r) 2 
Introducing new functions 

C{T) = V C 2A : e 3m * /4 (m(r , ) 1/2 + 2T"(T')- 1/2 ) , 

B(X) = (X") 2 (X')- 5/2 ~ 2X ,,, (X , )- 3/2 , 
D(X) = 2XB + 8X"(X')- 1/2 , 

Q(T) = ^p^e imt ^{m{T r ) l ' 2 T + 2(T'y x l 2 T"T - 4(T') 3/2 ), 

R(X) = BX 2 + 8X ,, (X , ) _1/2 X - -8(X') 3/2 , 
one can write the last equation as 

C(T)X + D(X)T + B{X)T 2 + Q(T) + R(X) = 0. (34) 

It is possible consider (34) as condition of orthogonality of two vector functions 
Z=(C,T,T 2 ,Q,1), W = (X,D,B,1,R). 
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Denote by p(Z) and p(W) the number of linearity independent functions 
among C, T, T 2 , Q, 1 and X, D, B, R, 1 respectively. From orthogonality condi- 
tion it follows that p(Z) + piW) < 5. It is possible to show that if T' ^ then 
p(Z) = 3 and piW) = 2. In this case we have 

D(X) = a x X + 6i, B(X) = a 2 X + 6 2 , = a 3 X + 6 3 , 

with a,i, bi £ It*. Because of (34) and definition of the functions D,B,Rwe obtain 
equations 

(X') 3 = (c 3 X 3 + c 2 X 2 + Cl X + c ) 2 , (35) 

(T'f = A(-c,T 3 + c 2 T 2 - Cl X + c ) 2 , (36) 

where C3,C2,ci and Co are arbitrary constants, A = (— 2r) 1 / 3 . 

The solutions of (35) and (36) can be expressed in the terms of the Weier- 
strass function p [ |I5j| . Indeed, one can write (35) and (36) as 

(X') 3 = (c 3 (X - ai )(X - a 2 )(X - a 3 )) 2 , (38) 

(Tf = A(-c 3 (T + ai)(T + a 2 )(T + a 3 )) 2 . 
Replacing X — a\ by 1/Y in (38) yields 

(Y l f + B 2 (Y-b 1 ) 2 (Y-b 2 ) 2 = J 

where B = c 3 («2 — «i)(a 3 — «i), &i = — - — , 62 = — z — • Introducing new function 



Z such that Z 3 = B(Y — b\)(Y — b 2 ), we obtain equation 

(Zy + ^ + fi \tf .Q. (39 ) 

The solutions of the last equation are expressed in the terms of the Weierstrass 
function p. Applying the above process to (36), we obtain equation such as 

(39). 

We shall omit here other cases and discuss briefly in the next section two- 
dimensional equation. 



4. Two-dimensional equation. 

We consider here the fast diffusion equation 



w f = Aln(u). (40) 
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Some applications of this equation can be found in JT6] . Galaktionov [[HJ used 



invariant subspaces to find some solutions of (40). If we set u = 1/v, we obtain 

v t = v 2 A\n(v). (41) 

It is easy to check that one of elementary solutions of (41) is travelling wave 
given by 

v = 1 + cexp(mx + ny — (m 2 + n 2 )t), (42) 

where c, m and n are arbitrary constants. Obviously, this is invariant solution. 
On the other hand, this solution satisfies differential constraints 

v x — mv + m = 0, 

v y — nv + n = 0. 

It is possible to find other differential constraints that are linear with respect 
x, y and v. 

It can be shown that the differential constraints 

ci — cotanit) ( 2 2\\ 

v x H 2~—^ — [v-xcq- yci + t(c + cj = c , 

c ~r c l V / 

c + dtan{t) ( 2 2 A 

% 2^ 2 I « - SCO - 2/Ci + t(C + Cj J = Ci 

are compatible with the equation (41). Here cq and c\ are arbitrary constants. 
The solution of (41) corresponding to these constraints is 

v = cqx + c\y — (cq + c 2 )t + C2 cos(t) exp(mix + m2?/ + m^t). 

Here co,ci,C2 are arbitrary constants and 

c t - Ci C\t + c / , \ 

mi = —2 2", 77T-2 = — 2 2"' 7713 = —tan(t). 

Cq + Ci Cq + C X 



We can derive other explicit solutions using invariant subspaces [JTt| or linear 



differential constraints. For example, from [[TTj one may extract the following 
representation 

v = so + s\ cos(x) + S2 sin(x) + S3 exp(y) + s^exp(—y), 

where functions s^-(t) satisfy ordinary differential equations 

s + s\ + s\- 4s 3 s 4 = 0, (43) 
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4 + sis = 0, 




4 + s 2 s = 0, 


(44) 


4 - S 3 S = o, 


(45) 


4 ~~ S 4S0 = 0. 





Because of (44) and (45) we find 

4/4 + S 3/S2 = o. 

This yields 

s 2 = c 2 /s 3 , c 2 G R. 
By arguments similar to that used above we have 

si = ci/ s 3 , s 4 = c 4 s 3 , ci, c 4 G i?. 

Substituting this into (43) leads to 

4 + ( C l + C 2) S 3 2 ~~ 4C4S3 = 0. 

From (45) we express the function so and obtain 

(ln S3 )" = ^-& S 3 2 , (46) 

with a = 4c4, b = cf + c|. 

Setting a = b = 1 one can derive two elementary solutions 

S3 = tanh(i), S3 = tan(t) 

In general, the solutions of (46) can be expressed in terms of elliptic functions. 
It is easy to obtain the correspondent function u. 
The more complicated representation is 

v = so + si cos(2x) + s 2 sin(2x) + S3 exp(2y) + S4 exp(— 2y)+ 

S5 sin(x) exp(y) + sq sin(x) exp(-y) + S7 cos(x) exp(y) + s$ cos(x) exp(— y), 

where Si are functions which satisfy some ordinary differential equations. The 
special case of this representation was found in |TT][ . 

It is important to note that the equation (40) is invariant under infinite- 
dimensional algebra of symmetry 0. Some solutions of (40) were obtained 
by means of these symmetries in |J. We shall describe other method of using 
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symmetry. It is convenient to apply the complex conjugate variables z = x+iy, 
z = x — iy. Thus, we can write the equation (40) as 

1 d * u 

Ui = lOzOz- (47) 
It is easy to check that (47) is invariant under the transformation 

z' = A{z), z' = B(z), u' = u/(A z B g ), 

where A(z) and B{z) are arbitrary functions. In other words, if the function 
f(t,z,z) is a solution of (47) then f(t, A(z), B(z))A z Bz also satisfies (47). 

For example, if we set m = n = 1 then from (42) we can construct the 
solution of the equation (40) 

A z Az 

u = 



1 + cexp(^ + A - 2t) ' 
where A is an arbitrary function of z and A is the complex conjugate function. 



5. Conclusions. 

In sections 2 and 3 we have shown how the method of the linear determining 
equations can be applied to find explicit solutions to nonlinear diffusion equa- 
tions. We have found exact solutions of these equations, using only the simplest 
solutions of the linear determining equations. It is interesting to find solutions 
of the linear determining equations depending on derivatives of higher orders. 
A.Shmidt |T7|, [T^] applied this method to another parabolic equations and some 
systems; application to the elliptic equation is discussed in 

In section 4 we have considered the two-dimensional equation. Applications 
of systems of the linear determining equations to multi-dimensional equations 
briefly discussed in Using results of section 3 one can find the following 



representation 

u = {a + be mt ' 2 ) 2 . (48) 

of solution of the equation 

u t = A(w 1 / 2 ) + mu + nu 1 ^ 2 , m,n G R, 
where the functions a(x, y) and b(x, y) must satisfy the system 

Aa = ma 2 + na, 
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Ab = mab + nb. 



It is easy to show that the differential constraint 

u tt = u t /2u J r mu t /2 

leads to the representation (48). The interesting reductions of some diffusion 
equations in several independent variables can be found in [|2D| , [2T| . It is impor- 
tant to explain these reductions by means of differential constraints. 
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